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REGULARITY CONDITIONS IN THE REALISABILITY 
PROBLEM WITH APPLICATIONS TO POINT 
PROCESSES AND RANDOM CLOSED SETS 
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Universite Lille 1, Universite Paris Descartes 
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By Ilya Molchanov''" 
University of Bern 

We study existence of random elements with partially specified 
distributions. The technique relies on the existence of a positive ex- 
tension for linear functionals accompanied by additional conditions 
that ensure the regularity of the extension needed for interpreting it 
as a probability measure. It is shown in which case the extension can 
be chosen to possess some invariance properties. 

The results are applied to the existence of point processes with 
given correlation measure and random closed sets with given two- 
point covering function or contact distribution function. It is shown 
that the regularity condition can be efficiently checked in many cases 
in order to ensure that the obtained point processes are indeed locally 
finite and random sets have closed realisations. 



1. Introduction. Defining the distribution of a random element in 
a topological space X is equivalent to specialising the expected values for 
all bounded continuous functions g(^). These expected values define a linear 
functional $(g) = Eg(^) on the space of bounded continuous functions g : 
A' I— 7- M. It is well known that a functional $ indeed corresponds to a random 
element if and only if $ is positive (i.e. $(g) > if g is non-negative) and 
upper semicontinuous (i.e. $(gn) 4 if gn i 0), see e.g. [35]. 

Below we consider the case of functional $ defined only on some functions 
on X and address the realisability of i.e. the mere existence of a random 
element ^ such that $(g) = Eg(^) for g from the chosen family G of functions. 
The uniqueness is not on the agenda, since typically the family G will not 



*The first author is grateful to the University of Bern for hospitality. 

^Supported by Swiss National Foundation Grant No. 200021-126503. The revision was 
carried over while the second author held the Chair of Excellence at the University Carlos 
III of Madrid supported by the Santander Bank. 

AMS 2000 subject classifications: Primary 60D05; secondary 28C05 47B65 60G55 

Keywords and phrases: point process, correlation measure, random closed set, two-point 
covering probability, contact distribution function, realisability 

1 



2 



R. LACHIEZE-REY AND I. MOLCHANOV 



suffice to uniquely specify the distribution of ^. A classical example of this 
setting is the existence of a probability distribution with given marginals, see 
[8]. The present paper focuses on some geometric instances of the problem. 
We will see that in most cases the answer to the existence problem consists 
of the two main steps. 

1. (Positivity) Checking the positivity condition on <I> — in most cases 
this requires checking a system of inequalities, which is a serious (but 
unavoidable) computational burden. 

2. (Regularity) Ensuring that the extended functional is regular (namely, 
upper semicontinuous) and so defines a u-additive measure. 

The first step ensures that it is possible to extend functional $ positively 
from a certain family of functions to a wider family. In this work we put 
the emphasis on the latter step — checking the regularity condition, leaving 
aside the computational difficulties arising from validating the positivity 
assumption. 

The use of positive extension techniques (that go back to L.V. Kan- 
torovitch) in the framework of stochastic geometry was pioneered by T. Kuna, 
J. Lebowitz and E.R. Speer [12] in application to point processes, which 
greatly inspired the current work. In this paper we establish the general 
nature of an idea proposed in [12] and show how it leads to various further 
realisability results. The new idea is to introduce an additional function, 
what we call the regularity modulus, and to formulate sufficient and nec- 
essary conditions in terms of a positive extension of a functional onto the 
linear space containing the regularity modulus and requiring only a priori 
integrability of the regularity modulus. 

We concentrate on two basic examples of the realisability problem: the 
existence of point processes with given correlation (factorial moment) mea- 
sure and the existence of a random closed set with given two-point cov- 
erage probabilities or contact distribution functions. The introduction to 
the realisability issue for point processes is available in several papers by 
T. Kuna, J. Lebowitz and E.R. Speer [11, 12], see also Section 3 of this 
paper. The realisability problem for random closed sets has been widely 
studied in physics and material science literature, see [7, 15, 30, 32, 33] and 
in particular the comprehensive monograph by S. Torquato [31] and a re- 
cent survey by J. Quintanilla [22]. If ^ is a random closed set (see Sec. 4 for 
formal definitions) in a locally compact Hausdorff second countable space 
X, its one-point covering functions is defined by 
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It is easy to characterise all one-point covering functions of random closed 
sets as follows. 

Theorem 1.1. A function px, x € X, with values in [0,1] is the one- 
point covering function of a random closed set if and only if p is upper 
semicontinuous. 

The upper-semi continuity of the one-point covering function of a random 
closed set is a straightforward consequence of the upper semicontinuity 
property of the capacity functional of a random closed set, see [20, Sec. 1.1.2]. 
Conversely, the function p from the theorem is realised (for instance) as the 
one-point covering function of the random set ^ = {x : Px ^ v} where v is 
a uniformly distributed variable (the details are left to the reader). 

It is considerably more complicated to characterise two-point covering 
functions 

Px,y = 'P{x,y G C} , x,yeX. 

In view of applications to modelling of random media it is often assumed that 
^ is a stationary set in MJ^, so that the one-point covering function is constant 
and the two-point covering function px^y depends only on x — y. Since a 
random closed set can be considered as an upper semicontinuous indicator 
function, the realisability problem for the two-point covering function can 
be rephrased as follows. 

Characterise covariance functions of (stationary) upper semicontinuous ran- 
dom functions with values in {0, 1}. 

These covariances are obviously a subfamily of positive semi-definite func- 
tions. Without the upper semicontinuity requirement, this problem, of com- 
binatorial nature, was solved by B. McMillan [19] and L. Shepp [26, 27] using 
the extension argument from [8]. More exactly, they normalised indicators 
by letting them take values -|-1 or —1 and assumed that the mean is zero. 
Their result does not rely on the topological structure of the underlying 
space and so does not necessarily lead to an upper semicontinuous indicator 
function. 

Example 1.2. Let px^y = \ and let Px = \ for ah x,y G M. While 
this two-point covering function corresponds, e.g., to the indicator field with 
independent values, it cannot be obtained as the two-point covering function 
of a random closed set, see Proposition 4.4. 

Even leaving aside the upper semicontinuity property, the McMillan- 
Shepp condition involves a family of corner-positive matrices, which is poorly 
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understood. As a result, its practical use to check the realisability for ran- 
dom media is rather limited. A number of authors have attempted to come 
up with simpler (but only necessary) conditions, see, e.g., [7, 17, 22, 32]. 
Another set of conditions for joint distributions of binary random variables 
is formulated in [25] in terms of the corresponding copulas. 

The realisability problem can be also posed for point processes in terms of 
their moment measures. In case of moment measures of arbitrary order it has 
been solved by A. Lenard [13, 14]. The case of moment measures up to the 
second order has been studied by T. Kuna, J. Lebowitz and E.R. Speer [11], 
whose recent paper [12] contains (among other results) a complete solution of 
this realisability problem for point processes with finite third-order moments 
and hard-core type conditions with fixed exclusion distance. The results of 
[12] can be extended to higher order moment measures, as was explicitly in- 
dicated there. Again, the positivity condition of [12] is extremely difficult to 
verify, even more complicated than the original condition for point processes 
because of new polynomial functionals involved in the positivity condition. 

The paper is organised as follows. Section 2 presents a series of general 
results on regular extensions and also invariant extensions (relevant for the 
existence of stationary random elements). These results form the theoretical 
backbone of our study, and are new even in the abstract setting of extending 
general positive linear functionals. 

Section 3 presents a number of realisability conditions for correlation mea- 
sures of point processes that considerably extend the results of [12] by re- 
laxing the moment and hardcore conditions. One of our most important 
results is Theorem 3.3 that shows how to split the positivity and regularity 
conditions, so that the latter can be efficiently checked. The importance of 
the packing number in relation to realisability conditions for hard-core point 
processes is also explained. 

Section 4 deals with the realisability problem for two-point covering prob- 
abilities of random sets. The closedness of the corresponding random set 
can be ensured by imposing appropriate regularity conditions. Section 5 ad- 
dresses a further variant of the realisability problem that involves contact 
distribution functions of random sets. 

The notational convention is that the carrier space is denoted as X (e.g. 
M,M'^), points in the carrier space are x,y, subsets of carrier spaces are 
denoted by capitals X, Y, F (while Y is reserved for counting measures iden- 
tified with corresponding support sets), the families of sets (or families of 
counting measures) as X,Af,T (while in Section 2 X denotes also a rather 
general space), and random element in these spaces (random sets or point 
processes) as ^, real functions acting on J- are g, v and families of such 
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functions are G, E, V, a functional on G, E,V is denoted by <I>, real numbers 
are denoted by t, r, A, while c denotes a generic constant and at the same 
time the corresponding constant function. 

2. Extending positive functionals. Fundamental results about the 
extension of positive operators form the heart of our main results, and are 
necessary to understand the machinery of the proofs. Nevertheless, the re- 
sults of the subsequent sections can be understood without Section 2, with 
the exception of Definition 2.5. 

2.1. General extension theorems. Consider a vector lattice E, i.e. a linear 
space with a partial order and such that for any vi , V2 G E their maximum 
vi V V2 also belongs to E. The absolute value |v| of v is defined as the sum 
of V V and (-v) V 0. 

Let G be a vector subspace of E, which is not necessarily a lattice itself, 
i.e. G may be not closed with respect to the maximum operation. We say 
that G majorises E if each v G E satisfies |v| < g for some g € G. A real- 
valued functional $ defined on E (resp. G) is said to be positive if <I>(v) > 
whenever v > and v € E (resp. v G G). A functional defined on E is said to 
be an extension of <^ : G i— t- M if it coincides with <1> on G. The extended $ 
is always denoted by the same letter. The following result about extension 
of positive functionals goes back to L.V. Kantorovich. 

Theorem 2.1 (see [1], Th. 8.12 and [34], Th. X.3.1). Assume that G is 
a majorising vector subspace of a vector lattice E. Then each positive linear 
functional on G admits a positive extension on the whole E. 

If G is a lattice itself, then it is possible to gain much more control over the 
extension of $, e.g. a continuous functional admits a continuous extension, 
see [34, Sec. X.5]. On the contrary, very little is known about regularity 
properties of the extension if G is not a lattice. 

In the following we assume that G and E are families of functions g on a 
certain space X. If G contains constant functions, the positivity of $ over G 
can be equivalently formulated as 



xex 

This equivalence is a particular case of the following result for x = (replace 
g with -g in (2.2)). 



(2.1) 




Proposition 2.2. Assume that vector space G contains constant func- 
tions and denote by G\M. the family of non-constant functions from G. If x 
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is any non-negative function on X , then a linear functional ^ on Q admits 
a positive extension on G + Mx with <I?(x) = r if and only if 

(2.2) r = sup $(g) = sup inf [x{X) - g{X)] + ^>(g) < oo . 

Proof. Since every element of G can be written c + g with g G G \ M and 
c G M, the left-hand side of (2.2) equals 

r = sup sup c + <^(g) = sup Cg + <l>(g) , 

g6G\Rc6M: c+g<x gGG 

where Cg = mix^xix ~ ?>){X) is the largest c such that c + g < x-, which 
yields the equality in (2.2). 

The necessity of (2.2) is straightforward because r < $(x) < oo. For 
the sufficiency, assume that (2.2) holds. The proof consists in checking that 
assigning the value ^{x) = yields a positive extension on G + Mx- Let us 
first prove that <^ is positive on G. If some g < satisfies <I>(g) > 0, then 
$(tg) t CO as t — )• oo whereas tg < x, which contradicts (2.2). 

Let g + Ax > for A / and g G G. If A > 0, then — A^^g < x, whence 
$(-A-^g) < r and ^'(g + Ax) > -Ar + A^>(x) = 0. If A < 0, -A'^g > x 
whence — A~^g is larger than any g' < x, and 

$(-A"^g) > sup <i>(g') = r 

g'eG,g'<x 

by monotonicity of <I> on G. Hence ^'(g + Ax) > — Ar + A<I>(x) = 0. □ 

The advantage of the latter condition in (2.2) consists in the explicit 
reference to the space X where random elements lie instead of checking the 
inequality g < x- 

2.2. Regularity conditions and distributions of random elements. Let E 
be a certain family of functions v : ^ i— > M defined on a space X with lattice 
operation being the pointwise maximum and the corresponding partial order. 

Theorem 2.3 (Daniell, see Sec. 4.5 [3] and Th. 14.1 [10]). Let a vector 
lattice E consist of real-valued functions on X and let E contain constants. 
If ^ is a positive functional on E such that $(v„,) J, for each sequence 
v„ J, and <I>(1) = 1, then there exists a unique random element ^ in X , 
measurable with respect to the a-algebra generated by all functions from E, 
such that $(v) = Ev(^) for all v G E. 
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In view of the positivity of <I>, the condition imposed on <I> is equivalent 
to its upper semicontinuity on E. In this paper, we start with a functional <I> 
defined on a vector sub-space G C E and discuss the existence of a random 
element ^ e X such that <I)(g) = Eg(^) for all g G G. In this case ^ is said 
to be realisable as a probability distribution on X. 

Assumption 2.4. The vector space G of functions on X contains con- 
stants and, for each gi,g2 € G, there exists a g € G such that (gi V g2) < g. 

From now on assume that A' is a topological space. 

Definition 2.5. Given a vector space G of functions on X , a regularity 
modulus on A" is a lower semicontinuous function x : [0, oo] such that 

(2.3) 'H^ = {X^X: x{X) < g{X)} 

is relatively compact for each g G G (if all g G G are bounded, x is a 
regularity modulus if and only if it has compact level sets). 

Examples of regularity moduli are given in Sections 3 and 4. A function 
V : <Y I— >• M is said to be x-i"^9ular if v is continuous on Tig for each g in G. 
Each continuous function is trivially x-regular. The proof of the following 
central result is based on the ideas from the proof of [12, Th. 3.14]. It should 
be noted that our result entails not only the realisability, but also provides 
a bound for the expected value of the regularity modulus. 

Theorem 2.6. Consider a vector space G of functions on X satisfying 
Assumption 2.4 and such that each g from G is x-'f^Sular for a regularity 
modulus X- Let ^ he a linear functional on G with $(1) = 1. Then, for any 
given r > 0, there exists a Borel random element in X such that 



(2.4) 



if and only if 



fEg(e) = <I>(g) forallgeG, 
1ex(0 < r, 



(2.5) sup <I>(g) < r . 

gGG,g<x 

Proof. Condition (2.5) is necessary because g < X implies $(g) = 
Eg(0 < ExiO < r. 

SufRciency. Let E be the family of all x-regular functions v that satisfy 
V < g for some g G G. Each function v G E is Borel measurable. Note that E 
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contains all bounded continuous functions that generate the Baire fi-algebra 
on X being in general a sub-c-algebra of the Borel one. For each vi, V2 € E, 
the function vi V V2 is x-regular and is majorised by gi Vg2, where gi, g2 € G 
majorise vi and V2 respectively. In view of Assumption 2.4, E is a lattice. 

Without loss of generality assume that the supremum in (2.5) equals r. 
By Proposition 2.2, $ is positive on G and can be positively extended onto 
G + Mx with $(x) = T, and further on to E + Mx by Theorem 2.1. It remains 
to prove that the obtained extension satisfies conditions of Theorem 2.3. For 
that, we use an argument similar to that of [12]. First restrict the obtained 
functional $ onto E. Assume that x is strictly positive. Consider a sequence 
{v„, n > 1} C E such that v„ J, 0. For each n, let g„, be a function of G 
such that v„ < g„. Take e > 0. Then /C„ = {X : v„(X) > ex{X)} is a 
subset of relatively compact T-Lg^/e-, since x is a regularity modulus. Since 
v„ is continuous on Tig„/e, the set /C„ is closed and therefore compact. The 
pointwise convergence v„ J, yields that n„/C„ = (recall that x is strictly 
positive). Since {ICn} is a decreasing sequence of compact sets, /C„o = for 
some no, whence v„(X) < ex(^) for sufficiently large n. The positivity of 
$ on E + Mx implies $(v„) < e<&(x) = ^f, whence $(v.n) J, 0. Theorem 2.3 
yields the existence of a random element ^ in such that <I*(v) = Ev(^) for 
ah V G E. 

Since x is lower semicontinuous, it can be pointwisely approximated from 
below by a sequence {vn} of non- negative continuous functions. Then v„ = 
min(n, v„) belongs to E and also approximates x from below, so that Ev„(^) = 
*^(vn) < ^(x) = while the monotone convergence theorem yields 

Ex(C) = lim Ev„(C) <r. 

n^oo 

If X is not strictly positive, it suffices to apply the above argument to x' = 
1 + X and use the linearity of <I>. □ 

Condition (2.5), equivalent to (2.2), is expressed solely in terms of the 
values taken by $ on G, and therefore yields a self-contained solution of 
the realisability problem. It is not easy to check in general, but if x can be 
approximated by functions Xn G G, n > 1, then it is possible to "split" (2.5) 
into the positivity condition on $ and the uniform boundedness of $(xn)) 
n > 1. This idea is used successfully in several different frameworks, which 
justify the abstract setting of Theorem 2.6: in Section 3 for point processes 
(see Theorem 3.1), in Section 4.4 for random closed sets (see Theorem 4.9), 
and in [5] in the framework of random measurable sets with the regularity 
modulus being the perimeter of a set. 
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The realisability problem is particularly simple if X is compact and G con- 
sists of continuous functions. Then, for identically vanishing Xi Theorem 2.6 
yields the following result, which is similar to the Riesz-Markov theorem, 
see [10]. 

Corollary 2.7. Let X be a compact space with its Borel a-algehra. 
Consider a vector space G containing constants such that each g G G is 
continuous and a map $ : G iH> M such that $(1) = 1. Then there exists a 
random element £^ in X such that Eg(^) = $(g) for a// g € G if and only if 
^ is a linear positive functional on G. 

2.3. Passing to the limit. The following result shows that the family of 
all random elements that realise $ in the sense of (2.4) is weakly compact. 

Theorem 2.8. Assume that G satisfies Assumption 2.4 and consists of 
continuous functions on a Polish space X with regularity modulus x- Let <I> 
be a linear positive functional on G . Then the family lUt of all Borel random 
elements ^ that satisfy (2.4) for any given r > is compact in the weak 
topology. 

Proof. Since x is a regularity modulus, the set Hr/e is compact. By 
Markov's inequality, 

P{UK/e} = P{xiO>r/e}<e, 

for all C e so that M is tight. 

Let {£,n,n' > 1} be random elements from SOT. Assume that S^n converges 
weakly to some ^. Without loss of generality assume that the ^„'s are defined 
on the same probability space and converge almost surely to ^. Since x is 
non-negative, Fatou's lemma yields 

r > liminf Ex(^n) > Eliminf x(e«) > Ex(lim^„) = Bx{0 , 

where the lower semicontinuity of x also has been used. 

Take an arbitrary g € G and define Tixg as in (2.3). Let g'''(-^) = 
max(g(X),0) be the positive part of g. Then, for A > 0, 

Eg+(e„) = Eg+(e„) + Eg+(e„) . 

Since g is continuous, T-Lxg is closed (and compact), so that if ^„ € Tixg for 
infinitely many n, then also ^ € Hxg. Furthermore, Ag and also g itself, are 
continuous and bounded on Tixg, so that Fatou's lemma yields 

limsupEg+(e„) < Elimsup(g+(^„) ^eW^g) 

<Eg+(0 %«,^<Eg+(0. 
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Thus 

limsupEg+(^„) < + Eg+(0 < ^ + Eg+(e) . 

Since A is arbitrary, 

limsupEg+(en) <Eg+(0. 

Since g"*" is non-negative, Fatou's lemma yields that Eg"'"(^„) Eg+(^). By 
applying the same argument to the function (— g), limEg(^„) = Eg(^), so 
that Eg(^) = $(g) for all g € G. Therefore, ^ e m. □ 

The following result concerns realisability of pointwise limits of linear 
functionals. Special conditions of this type for correlation measures of point 
processes are given in [12, Sec. 3.4]. 

Theorem 2.9. Let n > 1} be a sequence of linear positive func- 

tionals on a space G that satisfies the assumptions of Theorem 2.8. Assume 
that 

(2.6) liminf sup $„(g)<oo. 

" geG, g<x 

If $r!.(g) — ^ ^(g) for all g G G, then $ is realisable as a random element ^ 
satisfying (2.4) and .such that ^ is the weak limit of random elements realising 
for a subsequence n^- 

Proof. By passing to a subsequence, it suffices to assume that (2.6) 
holds for the limit instead of the lower limit. Let be a random element 
that realises If r is larger than the limit of (2.6), then P{^n ^ ^r/el ^ 
so that {(,n} is a tight sequence. Without loss of generality assume that ^„ 
weakly converges to a random element ^. 

The pointwise convergence of yields that Eg(^„) — )■ <&(g) for all g € G. 
Now the arguments from the proof of Theorem 2.8 can be used to show that 
Eg(Cn) Eg(^), so that Eg(,^) = <&(g) for ah g € G, i.e. ^ indeed satisfies 
(2.4). □ 

2.4. Invariant extension. Consider an abelian group G of continuous 
transformations acting on X. For a function v on A", define 

{ey){x) = y{ex), eee,xex. 

A functional <l> is said to be 0-invariant if, for each € and v from the 
domain of definition of ^{9v) is defined and equal to $(v). 
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A Borel random element ,^ in A" is said to be Q -stationary if, for each 
6 G Q, 6^ has the same distribution as ^. A variant of the following result 
for correlation measures of point processes is given in [12, Th. 4.3]. 

Theorem 2.10. Assume that G is a Q-invariant space satisfying As- 
sumption 2.4 and consisting of x-regw/ar functions. Furthermore, assume 
that at least one of the following conditions holds: 

(i) G consists of continuous functions and x pointwisely approximated 

from below by a monotone sequence of functions g„ G G, n > 1. 

(ii) X 'is Q-invariant. 

Let ^ be a Q-invariant functional on G. Then, for every given r > 0, there 
exists a Q-stationary random element ^ in X satisfying (2.4) if and only if 
(2.5) holds. 

Proof. As in [12, Prop. 4.1], the proof consists in checking hypotheses of 
the Markov-Kakutani fixed point theorem. Let 9Jt be the family of random 
elements ^ that realise $ on G, and satisfy Ex(^C) < r for every € 0. The 
family 9Jt is easily seen to be convex with respect to addition of measures, 
it is compact by Theorem 2.8, and 0-invariant, since ^ is ©-invariant on G. 
It remains to prove that 9K is not empty. 

(i) In view of (2.5), it is possible to extend positively onto G + so 
that Ex(,^) < r. The 0-invariance of <I> on G together with the monotone 
convergence theorem imply that Ex(^0 = '^xiCj ^ '"i whence € 971. 

(ii) By Proposition 2.2, we can extend $ positively onto the 0-invariant 
vector space V = G+Mx- Since $ is 0-invariant on G, we have ^{9{g+tx)) = 
^{9g) + t(^{6x) = <^(g + tx) for g + tx ill V, whence $ is 0-invariant on V. 
According to [28, Th. 3], $ admits a positive ©-invariant extension to the 
space E -|- Mx, defined like in the proof of Theorem 2.6. The restriction of 
the obtained functional onto E corresponds to a random element ^'m. X that 
verifies (2.4) and satisfies E(6lv)(0 = ^(6'v) = ^>(v) = Ev(0, 6* G 6, for v 
in E. Since E contains all bounded continuous functions on X , 9^ and ^ are 
identically distributed for all ^ € 0. □ 

3. Correlation measures of point processes. 

3.1. Framework and main results. Let M be the family of locally finite 
counting measures on a locally compact separable metric space X. We denote 
the support of y G 7\A by the same letter Y , so that x £Y means y ({x}) > 1. 

Equip M with the vague topology, see [2] . A random element ^ in Af with 
the corresponding Borel cr-algebra is called a point process. Denote by A/q 
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the family of simple counting measures, i.e. those which do not attach mass 
2 or more to any given point. If ^ is simple, i.e. € A/q a.s., then ^ can be 
identified with a locally finite random set in X, which is also denoted by ^. 

For a real function /i on X x X and counting measure Y = - given 
by the sum of Dirac measures, define 

gh{y) = ^ h{xi,Xj), 

whenever the series absolutely converges, the empty sum being 0. Note that 
the sum in the right-hand side is taken over all pairs of distinct points from 
the support of Y, where multiple points appear several times according to 
their multiplicities. The value gh{Y) is necessarily finite if h is bounded and 
has a bounded support. The value gh(Y) is termed in [12] the quadratic 
polynomial of Y, while polynomials of order n > 1 are sums of functions of 
n points of the process, and are constants if n = 0. 

Let G be the vector space formed by constants and functions gh for h from 
the space of symmetric continuous functions with compact support. Note 
that G satisfies Assumption 2.4, since 

(Cl + ghi)y (C2 + ghi) < Ci V C2 + ghiV/i2 e G 

for all Cl, C2 S M and /ii, /i2 € "^o- Furthermore, each gh is continuous in the 
vague topology, and so is x-regular for any regularity modulus x- 

Assume that ^ has locally finite second moment, i.e. E^(A)^ is finite for 
each bounded A. The correlation measure p (also called the second factorial 
moment measure) of a point process ^ is a measure on X x X that satisfies 

(3.1) / h{x,y)pidxdy) =Egh{0 

JXxX 

for each h € see [2, Sec. 5.4] and [29, Sec. 4.3]. The left-hand side defines 
a linear functional $(g/i) on g/^ E G. 

Let X he & subset of M, which may be A/" itself. Given a symmetric locally 
finite measure p on X x X, the realisability problem amounts to the existence 
of a point process ^ with realisations from X and with correlation measure 
p, so that ^{gh) = EghiO for ah hG'tfo- 

By (2.1), the positivity of $ means 

(3.2) cl>(g,) > inf g,(y) 

for all /i € '^o- Then it is clear that the positivity of ^ is necessary for 
its realisability. If X is compact in the vague topology, then Corollary 2.7 
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applies and the positivity condition (3.2) is necessary and sufficient for the 
realisabihty of p. 

However, in general the positivity condition alone is not sufficient for the 
realisabihty, see [11, Ex. 3.12]. In the following we find another condition 
that is not directly related to the positivity, but, together with the positivity, 
is necessary and sufficient for the realisabihty. 

As an introduction, let us present our results for X being a subset of the 
Euclidean space W^. For e > 0, define 

xe(X)= Yl ll^-yll"''"'' ^e-^' 

x,y£Y, Xy^y 

which is later acknowledged as being a regularity modulus (see Definition 
2.5) if e 7^ 0. Note that XeO^) is infinite if Y has multiple points. The tools 
developed in this paper enable us to resolve the original realisabihty problem 
with a supplementary regularity condition involving Xe- 

Theorem 3.1. (i) Let be a compact subset of M.'^ without isolated 
points. A symmetric finite measure p{dxdy) on X x X is the correlation 
measure of a simple point process ^ C X such that Exo < oo if and only 
if ^ given by the left-hand side of (3.1) is positive and 

/ \\x — y\\~'^ p{dxdy) < oo . 

(a) Let p be a symmetric locally finite measure on x such that p{[A-\- 
x) X (B-\-x)) = p{A X B) for all x € M'^ and measurable sets A and B. Then 
there exists a simple stationary point process ^ with correlation measure p, 
such that 

Exo(C n C) < oo 

for every compact C dW^, if and only if ^ defined by (3.1) is positive and 
(3.3) / ||x - y\\~'^p{dxdy) < oo 

J BxB 

for some open set B. 

Proof. The first statement follows from Theorem 3.5 using the fact that 
the packing number Pt(K) of X is bounded by ct~'^ for all sufficiently small 
t. For (ii), apply Theorem 3.9(ii) noticing that the imposed condition is 
equivalent to (3.24). □ 
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In the remainder of this section, one can find a quantification of this result 
(i.e. how the left hand member of (3.3) controls the value of Exo(^ H C)) 
as well as generalisations for general metric spaces. The main argument 
used is a splitting method based on Theorem 2.6, the details are made clear 
in the proof of Theorem 3.3. Note that the packing number of the metric 
space appears as a crucial quantity to uncouple in this way the realisability 
problem, see Lemma 3.2. 

3.2. Moment conditions. The family of all counting measures with 
total mass at most A; on a compact space X is compact. Thus, a measure p 
on X X X is realisable as a point process with at most k points if (3.2) holds 
with X = X^. 

Assume that y is a finite counting measure. For a > 2 define 



The finiteness of Exq('^) amounts to the finiteness of the moment of order a 
for the total mass of ^. Since /i G is bounded by a constant c' and q > 2, 
the family 



{Y eM: Xo.{Y) < c + g;,(y)} C {y G AT : y(X)" < c + c'Y{Xf} 



consists of counting measures with total masses bounded by a certain con- 
stant and therefore is compact in the space Af. Hence Xa is a regularity 
modulus and so Theorem 2.6 yields the realisability condition 



of p by a point process ^ whose total number of points has finite moment of 
order a. Note that [12, Th. 3.14] provides a variant of this result assuming 
the existence of the third factorial moment of the cardinality of ^ (i.e. with 
a = 3) and for the joint realisability of the intensity and the correlation 
measures. The condition of [12, Th. 3.14] (reformulated for the correlation 
measure only) reads in our notation as c+^{gh)+br > whenever c+gh+bx3 
is non- negative on J\f. Noticing that b > 0, this is equivalent to the fact that 
c + $(g/i) < r whenever c + gh < X3) being exactly (3.4). If is a group of 
continuous transformations acting on X and p is G-invariant, then the point 
process ^ can be chosen 0-stationary by Theorem 2.10(ii). 
In order to handle possibly non-finite point processes ^ define 



xa{Y) = y(x) 



a 



Y eJ\f. 



(3.4) 



sup <I>(g) < oo 

gGG, g<Xa 
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for a lower semicontinuous strictly positive function /3 : X i— t- M and a > 2. 
By approximating (3 from below with compactly supported functions, it is 
easy to see that Xa,i3 is a regularity modulus. By Theorem 2.6 and (2.2), for 
any given r > 0, there is a point process ^ with correlation measure p such 
that Exq,/3(C) <r\i and only if p satisfies 



For a = 3, condition (3.5) is a reformulation of [12, Th. 3.17] meaning 
the positivity of $ on a family of positive polynomials that involve symmet- 
ric functions of the support points up to the third order. The realisability 
condition for 0-stationary random elements can be obtained by applying 
Theorem 2.10. 

3.3. Hardcore point processes on a compact space. Assume that X is a 
compact metric space with metric d. Let Me be the family of e-hard-core 
point sets in X (including the empty set), i.e. each Y € A/"e attaches unit 
masses to distinct points with pairwise distances at least e with a fixed e > 0. 
In this case no multiple points are allowed, i.e. Me C A/q. 

According to [6, 9], a subset X of simple counting measures TVo is relatively 
compact if and only if sup{Y (K) : Y € X} is finite and the infimum over 
Y £ X of the minimal distance between two points in Y n K is strictly 
positive for each compact set i^' C X. The hard-core condition yields that 
the number of points in any compact set is uniformly bounded, and so A4 is 
indeed compact. By Corollary 2.7, p is realisable as the correlation measure 
of an e-hard-core point process with given e > if and only if 



for all h G This result is formulated in [12, Th. 3.4], which essentially 
reduces to the positivity of $ over the family c + gh (in our setting) . 

In this paper we assume that the hardcore distance is not predetermined 
and the point process takes realisations from Ue^oMe, which coincides with 
A/q in case of compact X. Note that (3.6) is stronger than the positivity of 
$ on functions g/j defined on the whole family A/q and formulated as 



(3.5) 




(3.6) 



(3.7) 



If X does not have isolated points, then the infimum in (3.7) can be taken 
over J\f. This is seen by approximating a multiple atom with a sequence 
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of simple counting measures supported by points converging to the atom's 
location. 

In the following we use the (hard-core) regularity modulus of the form 

where ip : (0, oo) i-^- [0, oo] is a monotone decreasing right-continuous func- 
tion, such that 'ip{t) oo as t ^ 0. The compactness of X and the lower 
semicontinuity of ip imply that is lower semicontinuous on TVq. As shown 
below is a regularity modulus if ip grows sufficiently fast at zero. 

Let Pt(K) be the packing number of X, i.e. the maximum number of points 
in X with pairwise distances exceeding t, see [18, p. 78]. It is convenient to 
define the packing number at t = as i-b(X) = oo if X is infinite and 
otherwise let Po(^) be the cardinality of X. 

Lemma 3.2. Function x^ is a regularity modulus on Mq if 

(3.8) V(i)/-ft(X) ^ oo as tiO. 

Proof. In view of the compactness of X, it is possible to bound /i G "^o 
by a constant A, so that x!^ is a regularity modulus if 

nx = {Y(^N,: x^'^iY) < XY{Xf} 

is compact in A/q for each A > 0. For this, it suffices to show that the total 
mass of all Y G Hx is bounded by a fixed number and Hx C Afe for some 
e > 0. 

Let 7t(n) be the minimal number of pairs {xi,Xj) with i ^ j, such that 
Xi,Xj G Y and d{xi,Xj) < t over all counting measures Y of total mass n. 
Take t such that ^|;{t)/Pt(K) > A. If y(X) > n, then 

xi'iY) > Yl ^(*) 2d(x„x,)<* > lt{n)m . 

Xi,Xj£Y, i^j 

Therefore, 

-HxCiY: n~\{n)i;{t) < A} 

consists of Y with total mass uniformly bounded by fixed number nx - Indeed, 
by Lemma A.l, 

lim n-\{n) > lim n'^n ( " 1 ) = Pt(^)~^ ■ 

rt— >-oo n— >-oo \r^i[2L) J 



REALIS ABILITY PROBLEM 17 

Choose e > so that ip{t) > Xnj^ for t < e. For Y £ Hx and any x,, Xj € Y , 
^{d{xi,Xj))<x':p{Y)<\nl, 
whence d(a;j, Xj) > e. Thus T-Lx C Me-, so "Ha is relatively compact. □ 

The following theorem shows that the realisability condition can be split 
into the positivity condition (3.7) on the linear functional $ and the regu- 
larity condition (3.9) on the correlation measure, so that the latter can be 
easily checked. Such a split is possible because the regularity modulus 
can be approximated by functions from G. 

Theorem 3.3. A locally finite measure p on % x ^ is the correlation 
measure of a simple point process ^ such that ^x^iO — ^ /^'^ some r > 
with tp satisfying (3.8) if and only if (3.7) holds and 



(3.9) / ilj{d{x,y))p{dxdy) <r. 

Proof. Necessity. The definition of the correlation measure implies that 
i^{d{x,y))p{dxdy) =Bx':p'iO <r ■ 

X 

SufEciency. First assume that tp only takes finite values. The proof consists 
in checking (2.2), which is equivalent to (2.5). 
For each family of positive numbers {tg,g G G}, 

(3.10) sup inf [x{Y) - g{Y)] + <D(g) < sup inf [x{Y) - g{Y)] + $(g) . 

geG ^6M) gGG ^GMg 

The crucial step of the proof consists in the careful choice of tg > 0. 

Fix g € G. For t > 0, define V't(s) = ip{max{t, s)), s > 0. Since any Y G Aft 
does not contain any two points at distance less than t, x(^) = Siiai-^)- 
Therefore, 

(3.11) yigj^J^(^) - = ~ ^^^^^ • 
Our aim is to prove that 

(3.12) ^inf.^(g^, - g)(y) = ^n|^(gV.. " g)(^) , 

because then, since g^^ G G, the positivity of $ on G yields that (3.11) is 
not greater than ^'(g^i — g). Thus, (3.10) is bounded above by 

sup $(g^, -g) + ^{g)< sup $(g^, ) = / ip{d{x, y))p{dxdy) 
gGG t JXxX 
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by the monotone convergence theorem. 

The proof of (3.12) rehes on the proper choice for t (depending on g). 
Assume without loss of generahty g = g/j for /i G 'rfo with absolute value 
bounded by A > 0. By (3.8), there exists to such that V'(^o)/-fto(X) > A + 1. 
By Lemma A.l, there is uq such that for all Y with mass n > uq, the number 
of pairs of points of Y at distance at most Iq satisfies 



Choose t < to so that 'tp{t) > Xn^ and consider any Y € Mq \ Aft- If Y(K) < 
no, then 

while if y(X) > no, then 

g^,{Y) - g,(y) > g^,^(y) - g,(y) > nV(to)7to(>^) - An^ > o. 

Thus for Y ^ Aft, we have gip^iY) — ghiY) > 0. Therefore, the infimum of 
SVt ~&h, which is non-positive because zero is obtained for y = 0, is reached 
on Aft, and (3.12) is proved. 

Now assume that ■i{j{t) is infinite for t € [0,6) and finite on {S,oo) with 
5 > 0. If ip{t) oo as t -I S, then the above arguments apply with tQ > d 
chosen such that ip{to)/ P5(^) > A. 

Assume that ipi^) is finite. Let tpo{t) be a function satisfying (3.8) and 
finite for ah t > 0, e.g. ^o(i) = t"^-Pt(X). Define V*(0 = V'(i) for t > 5 
and let 'ip*{t) = tpo{t) + a for t G (fi,6) with a sufficiently large a, so that 
ifi* is monotone right-continuous, and xjj,* is a regularity modulus. Applying 
the previous arguments to ip* yields that there exists a point process ^ such 
that Ex5j,*(0 < r. Since r < oo, p vanishes on {{x,y) : d{x,y) < 5}, and so 
Ex^^,(0=Ex^^(e)<r. □ 

The following result is obtained by letting ■0 be infinite on [0,e) and 
otherwise setting it to zero. 

Corollary 3.4. A measure p on X x X zs the correlation measure of 
a point process ^ with ^ S Af^ a.s. if and only if (3.7) holds and p{{{x,y) : 
d{x,y) <e}) = 0. 



The following result yields a direct realisability condition for p without 
mentioning a regularity modulus. 
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Theorem 3.5. Let p be a locally finite measure on X x X, and fix any 
r > 0. Then there exists, for every r' > r, a simple point process ^ with 
correlation measure p, such that 

(3.13) E Yl ^'d(.,,x,)(X)<r', 
if and only if ( 3. 7) holds and 

(3.14) / Pd(x,y)(X)p(dxdy) <r. 

Proof. Necessity. Call /ij(x,y) = min(t, j^)(X)) for a; 7^ y € X and 
t > 0. Assume that ,^ realises p and satisfies (3.13). The monotone conver- 
gence theorem yields that 

/ Pd{x,j;) {^)p{dxdy) = lim Eg^, (0 < r' 

for every r' > r, whence (3.14) holds. 
SufEciency. Define a measure on M_|_ by 

p'{[a, h)) = p({(x, y) G X X X : a < d(x, y) < h}) . 

Fubini's theorem yields that 

PtiX)p'idt). 



Let {tk, k > 1} he a strictly decreasing sequence of numbers such that 

[ Pt{X)p{dt) < 2-^ 

J[0,tk) 

For m > 1, the function 



kPt{X) iftk+i<t<tk<t^,k>l, 
Pt{X) ift>t^ 



is monotone right-continuous and satisfies 'iprn.{t) / PtiX) ^ 00 as t — )> 0. Then 
'iljm{^{x,y))p{dxdy) = j ip„i{t) p' (dt) 

< I Pt(x)p'{dt) + ^'^'^ ^ ^ + ^'^'^ ■ 



k>m k>m 
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By Theorem 3.3, choosing m sufficiently large yields the realisability of p by 
a point process ^ satisfying 

Xi,XjG^, i^j k>m 

□ 

Remark 3.6. Let be a group of continuous transformations on X 
that leave p invariant, i.e. p{9A x 6B) = p{A x B) for all € and Borel 
A,B. Since the regularity modulus can be approximated from below 
by a sequence of functions from G, Theorem 2.10(i) is applicable and so 
the corresponding point process ^ in Theorems 3.3, 3.5 and Corollary 3.4 
can be chosen ©-stationary. If Q consists of isometric transformations, then 
Theorem 2.10(ii) is also applicable. 

3.4. Non-compact case and stationarity. Assume that X = R'^ and d(x, y) = 
\\x — y\\ is the Euclidean metric. Let tp he a positive right-continuous mono- 
tone function on such that ilj{t)t'^ — )■ oo as t ^ 0. Denote by Bn the 
open ball of radius n centred at 0. Given a known bound for the packing 
number in the Euclidean space [18, p. 78], Lemma 3.2 implies that is a 
regularity modulus on every n > 1. Define 

(3.15) X%{y)= P{x^,Xj)iJ{\\xi-X,\\) 

for a bounded lower semicontinuous strictly positive on X function (5. 

Theorem 3.7. Let p he a locally finite measure on x W^. 

(i) The measure p is realisable as the correlation measure of a point process 
^ that satisfies Ex^'ipiO — ^/ ^'^^ '^'^^^ ^/ (^■'^) holds and 

(3.16) / (3{x,y)'ilj{\\x -y\\)p{dxdy) <r. 

(ii) Fix r >0, let 

(3.17) rn= \\x - y\\~'^p{dxdy), n > 1, 

J Bn X Bn 

and let {/?„, n > 1} be a sequence of non-increasing numbers converg- 
ing to 0. Then the following assertions are equivalent. 
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(a) (3.7) holds and 

(3.18) ^ /3n.(rn+i - rn) <r < oo, 

n>l 

in particular every r„, n > 1, is finite. 

(b) For every r' > r there exists ^ with correlation measure p and 
such that 

(3.19) ^(/3„ - /3„+i)E \\x^-Xj\\-'' <r' . 

n>l Xi,Xj&Bn,i^j 

Proof. SufRciency. (i) The function x}^^ is a regularity modulus on A/q, 
since 

n^h = {YeX: xf^Y) < c + gh{Y)}, c G M, /i € , 

is compact in TVq. This follows from Lemma 3.2, which yields the compact- 
ness of the restriction of Y from 7ic,h onto any compact set C. Indeed, this 
family of restricted counting measures coincides with the family of simple 
counting measures supported by C such that X^{Y) < c/m + gh/m{Y), 
where m > is a lower bound of f3{x, y) for x,y G C. 

In order to apply Theorem 2.6 with the regularity modulus (3.15) and in 
view of (2.5) it suffices to show that 

(3.20) inf \xf4Y) - gh{Y)] + [ h{x, y)p{dxdy) < r 

for all h G Assume that h is supported by a subset of x i?„ for 
some n > 1. Then (3.20) holds by the same reasoning as in the proof of 
Theorem 3.3 applied to the compact space Bn (One might first consider 
only Y C Bn, and then note that the infimum over all Y S A/q is necessarily 
smaller). By Theorem 2.6, (3.16) implies the existence of a point process ^ 
with correlation measure p that satisfies 'Eix^^{^) < r. 

(ii) Define Y„, = {Bn x Bn) \ (S„_i x S„_i), n > 1 (with Bq = 0). For 
every n > 1, define the measure 

Pn([a> b)) = p{{{x, y) € BnX Bn: a< \\x - y\\ < b}), 

and let 

i^=sup|t>0: / s-''p'n{ds)<2-H, k>l. 

J[0,t) I 
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Since p'^+i > for every n > 1, for every k,n > 1 we have tj!"*"^ < t^. Let 
{m„, n > 1} be a non-decreasing sequence of positive integers so that 

(3.21) Yl Pn{rn - r„_i + A:2-'=) < r' . 

n>l k>mn 

Now define 



Since m„ < rrin+i, ipn+i ^ V'n foi' every n > 1. Function ipn satisfies 
ipn{t)t'^ — 7> oo as t ^ 0, whence, for every n > 1, is a regularity moduhis 
on counting measures supported by and 

M\\x - y\\)pidxdy) < [ Mt)p'nidt) , 



where 

p';{[a, b)) = p{{{x, y) G Y„ : a < ||x - y\\ < b}) < p'^{[a, b)) . 

Then 

' M\\x-y\\)p{dxdy)< [ r'^p'^(dt)+ [ Mt)p'nidt). 



+ 



We have 

t~'^Pn{dt) = I \\x- y\\-'^ p{dxdy) = r„ - r„_ 



1 



(B„xB„)\{B„_ixB„_i) 



and 



whence 



Ktr, 



Mi)pnidt) < j Mi)pnidt) < J2 ^'^ 



k>mn 



(3.22) / i;ni\\x-y\\)p{dxdy) <rn-rn-i+ k2~K 

Define ijj{x,y) = ipn{\\x — y\\) for x,y £ Y.„. Since '(/'n+i < V'n and functions 
i^m n > 1, are lower semicontinuous, the function ^/^ is lower semicontinuous 
on X R'^. Define I3{x,y) = /3n on Y„. Since /3n, ^^ > 1, decrease, /3 is a 
lower semicontinuous function on x M'^. Since V'n ^ ^fc for every k < n, 
the restriction of onto sets y C Bn is larger than X/3^^„) whence 
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is a regularity modulus on A/q. By Theorem 2.6, is realised by a point 
process ^ satisfying 

Exf40 < [ /3(x, y)p{dxdy) = y2Pn [ M\\x - y\\)p{dxdy) . 

Since t"'^ < 'i/'nl^) for each n and t > 0, 

m 

n=l 
m 

> hm EV(/3„-/3„+i)x^=jeni?n)) 



m— >oo 

n=l 



> J](/3„-/3„+i)E Yl 

n>l ijtj,Xi,Xj&Bn 



Using successively (3.22) and (3.21) 

5^/3n / V'n(||x-y||)p(dxdy) < ^/3„(r„-r„_i+ ^2"'=) < r' , 

n>l "^^^ n>l fc>m„ 

we arrive at (3.19). 

Necessity. For (ii) remark first that for x ^ y (zW^ 

P{x,y)\\x - y\\-'^ = ^(/3n - /3n+i) Ix,j/gb„ Ik - 2/11"'^ , 

n>l 

where /3 is defined in the sufficiency part of the proof. The function I3tjj in (i) 
and the function (x,y) i->- fi{x ^y)\\x — y\\~'^ in (ii) are lower semi-continuous 
and can therefore be approximated from below by compactly supported con- 
tinuous functions. The rest follows from the monotone convergence theorem 
similarly to the proof of necessity in Theorem 3.5. □ 

Remark 3.8. Remark for point (ii) that if each r„, n > 1, is finite, 
there always exists a sequence {fin} of sufficiently small numbers such that 
the right-hand side of (3.18) is finite. 

If the distribution of point process ^ is invariant with respect to the group 
of translations of W^, then ^ is called stationary. Its correlation measure 
p is translation invariant, i.e. p{{A + x) x [B + x)) = p{A x B) for all x € M'^ 
and so 

(3.23) p{A xB) = \^ [ [ p{dy)dx , 

J A JRd 
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where A is the intensity of ^ and p is a measure on Mf^ cahed the reduced 
correlation measure, see [24, p. 76]. 

Theorem 3.9. Let p be a locally finite measure on M'^, let j3 be a bounded 
lower semicontinuous strictly positive function on satisfying 

^{y) = / (3{x,x + y)dx <oo, yeM.'^, 

and let ip be a monotone decreasing non-negative function such that t'^il){t) 
oo. 

(i) p is the reduced correlation measure of a stationary point process ^ that 
satisfies ExJ^'j^(C) ^ r if and only if (3.7) holds and 

[ mi^iMMdy) < r . 

(a) p is realisable as the reduced correlation measure of a stationary point 
process ^ that satisfies (3.19) for some sequence n > 1} if and 
only if 

(3.24) / < oo 

Jb 

for some open ball B containing the origin. If j^^i \\y\\~'^p{dy) is finite, 
it is possible to let /3„ = n~'^~^ , n>l, for any fixed (5 > 0. 

Proof. It suffices to use (3.23) to confirm the conditions imposed in 
Theorem 3.7, see also Remark 3.8. In order to show that ^ can be chosen 
stationary, note that x}^'^ can be pointwisely approximated from below by a 
monotone sequence of functions from G, so Theorem 2.10(i) applies. □ 

3.5. Joint realisability of the intensity and correlation. Recall that the 
intensity measure pi of a point process ^ is defined from 

E^/i(2;j) = j h{x)pi{dx) , h e '^0,1 , 

where ^0,1 is the family of continuous functions on X with compact sup- 
port. A pair {pi,p) of locally finite non-negative measures on X and X x X 
respectively is said to be jointly realisable if there exists a point process ^ 
with intensity measure pi and correlation measure p. 
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Let Gi be the vector space formed by constants and functions 

for hi € ^0,1 and h € It is easy to see that Assumption 2.4 is verified in 
this case. The pair (pi,p) yields a hnear functional 

(3.25) <^{gh,,h) = hi{x)pi{dx) + h{x,y)p{dxdy). 

Jx JXxX 

The realisability of <1> by a point process ^ means that ^{ghj^h) = ^ghi,hiC)- 
Functional ^ is positive on Gi if and only if 

(3.26) ^(ghi,h) > mf gh„hiY), hi € <^fo,i, h e ■ 

Similar arguments as before apply and yield the joint realisability con- 
ditions. Consider the special case of stationary processes in X = M"^ with 
the reduced correlation measure p (see (3.23)) and intensity pi{dx) = Xdx 
proportional to the Lebesgue measure. 

Theorem 3.10. Let X be a constant, and let p be a locally finite measure 
o/M'^. Then there is a stationary point process ^ with intensity pi(dx) = Xdx 
and reduced correlation measure p if ^ given by (3.25) satisfies (3.26) with 
X = A/q and 

z\\~'^p{dz) < 00 

for some open set B containing the origin. 

Proof. It suffices to note that gh^^h is dominated by cgh for a constant 
c and follow the proof of (ii) in Theorem 3.7. The condition on p follows 
from (3.17) and (3.23). □ 

4. Realisability of covering probabilities for random sets. The 

nature of realisability problem changes with the choice of the family of sub- 
sets of a carrier space X taken as possible values for a random set. We start 
with the case when a random set is allowed to be any subset of X, where re- 
alisability results are available under minimal conditions, while the obtained 
random set lacks properties and might even not be measurable. In the re- 
mainder of this section we treat the case of random closed sets, a classical 
setting in stochastic geometry. Some examples of possible regularity moduli 
are presented, along with the corresponding realisability results that resem- 
ble those of [11] in the point processes setting. The framework of random 



26 



R. LACHIEZE-REY AND I. MOLCHANOV 



measurable sets with finite perimeter (in the variational sense), developed 
in the forthcoming paper [5], provides a compromise between regularity of 
the random set and the explicitness of conditions. 

4.1. Random binary functions. Let X be the family of all subsets of 
X identified with their indicator functions. Endow X with the topology of 
pointwise convergence and the corresponding cr-algebra. Since X is compact, 
Corollary 2.7 yields the following result. 

Theorem 4.1. Let G be a vector space that consists of continuous func- 
tions on X and includes constants, and let ^ be a map from G to M. Then 
there exists a random indicator function ^, such that $(g) = Eg(^) for all 
g G G i/ and only if ^ is a linear positive functional on G and ^(1) = 1. 

The key issue in applying Theorem 4.1 is the choice of the space G. 

Example 4.2 (One-point covering function). Let G be generated by 
constants c and one-point indicator functions gx{F) = I^eF, F £ X, for 
X G X. The positivity of a linear functional <^ : G i— )■ M together with 
$(1) = 1 means that '^(g^) G [0, 1] for all x £ H. Thus, a function px = 
^{gx) is a one-point covering function P{x G (,} for a random set ^ if 
and only if takes values in [0, 1]. Compare with Theorem 1.1, where the 
extra upper semicontinuity condition ensures that the corresponding random 
binary function is upper semicontinuous and so ^ is a random closed set. 

Example 4.3 (Covariances of random sets). Consider vector space G 
generated by constants and functions gx^y{F) = 1x,y^F for x,y G X. The 
values of a linear functional $ on G are determined by Px,y = ^{gx,y), 
x,y By (2.1), $ is positive on G if and only if 

n n 

(4.1) Yl (^ijP^^,^j ^ jnf^ «U ^x,,x,&F 

ij=l ij=l 

for all n > 1 and all matrices (ajj)^-^]^. In particular, if then 
(4.1) implies the non-negative definiteness of Px,y, G X. Note that the 
one-point covering probabilities are specified if px,y are given. 

4.2. The closedness condition. A random closed set ^ in a locally com- 
pact Hausdorff second countable space X is a random element that takes val- 
ues in the family X = T of closed subsets of X equipped with the a-algebra 
(called the EfFros a-algebra) generated by families {F £ F : Fr\K ^ for 
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all compact sets K. The distribution of a random closed set ^, is uniquely 
determined by its capacity functional 

T{K) = P{^ n K / 0} 

for all K from the family of all compact sets in X, see [16] and [20, Th. 1.1.13]. 

Theorem 4.1 ensures only the existence of a binary stochastic process 
with given marginal distributions up to a certain order. However, it is not 
guaranteed that the constructed stochastic process has upper semicontinu- 
ous realisations, which should be the case if this process is the indicator of 
a random closed set in a topological space X. If the carrier space X is finite 
(more generally, has a discrete topology), then this problem is avoided, since 
each random set is closed. Furthermore, the closedness issue can be settled 
in the following special case of two-point probabilities in the product form 
(and can be generalised for multi-point covering probabilities). The follow- 
ing result implies, in particular, that the random indicator function from 
Example 1.2 does not correspond to a random closed set. It also illustrates 
regularity problems arising in realisability problems for random closed sets. 

Theorem 4.4. Assume that H is a separable space. A function 

_ jpxPy ifx^y, 

Px,y — \ 

[Px ijx = y 

is a two-point covering function of a random closed set if and only if Px, 
X £ is an upper semicontinuous function with values in [0, 1] such that 
each point x with px € (0,1) has an open neighbourhood U such that py > 
only for at most a countable number of y £ U and the sum of py for y £ U 
is finite. 

Proof. SufEciency. Note that the set L = {x : px = 1} is closed by 
the upper semicontinuity of px- The separability of X and the condition of 
theorem imply that the set M = {x : Px £ (0,1)} is at most countable. 
The sufficiency is obtained by a direct construction of a random subset Z 
of M that contains each point x with probability px independently of all 
other points. It remains to show that the random set ,^ = Z U L is closed. 
Consider x £ M and its neighbourhood from the condition of theorem. Since 
"^Py < oo, only a finite number of y belong to Z and so they do not converge 
to X. Thus, X with probability zero appears as a limit of other points from 
^ unless X £ L and so belongs to ^ almost surely. 

Necessity. The function px = F{x £ S^} is upper semicontinuous, since ^ 
is a random closed set. The product form of the two-point covering function 
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implies that the capacity functional on two-point set is given by 

T{{x, y}) =Px+Py- PxPy ■ 
The upper semicontinuity property of the capacity functional yields that 

limsupr({x,y}) < px , 

y^x 

while the monotonicity implies that T{{x,y}) px as y ^ x. Thus — 
Px) — >■ as y — >■ X for all x. Unless px = 1, we have py — )• 0. 

Assume that Px > ^ and px„ > 0, where — >■ x and x„ 7^ x with ^ px„ = 
00. A variant of the lemma of Borel-Cantelli for pairwise independent events 
(see [4, Lemma 6.2.5]) implies that almost surely infinitely many points Xn 
belong to ^, so that x G ^ a.s. by the closedness of ^ and so px = 1. Thus, 
the sum of Px^ for each sequence {x„} in a neighbourhood of x is finite. 
This rules out the existence of uncountably many y with > in any 
neighbourhood of x. Indeed, then {y : Py > 1/n} is finite for all n, and so 
the union of such sets is countable. □ 

It is known that J- is compact in the Fell topology that generates the 
Effros (T-algebra, see [20]. However, Corollary 2.7 is not applicable, since 
functions ILx.y^F, F T, generating the vector space G, do not generate the 
Effros cr-algebra on J^. 

It is known [20, Th. 1.2.6] that the cr-algebra generated by G on the family 
of regular closed sets (that coincide with closures of their interiors) coincides 
with the trace of the Effros cr-algebra on the family of regular closed sets. 
However, the family of regular closed sets is no longer compact in the Fell 
topology. Furthermore, indicator functions are not continuous in the Fell 
topology, so it is again not possible to appeal to Corollary 2.7 or explicitly 
check the upper semicontinuity condition required in Daniell's theorem. 

In order to describe a useful family G of functionals acting on sets, consider 
a (7-finite measure on X and define 

(4.2) gh{F) = [ h{x,y)u{dx)u{dy) 

J FxF 

for all measurable F C X and h from the family 'rfo of symmetric continuous 
functions with compact support in X x X. A function Px,y, x, y G X, generates 
a functional acting on gh as 



(4.3) 



^{gh) = / Px,yh{x,y)u{dx)u{dy) 



REALISABILITY PROBLEM 



29 



The function px^y is said to be weakly realisable as the two-point covering 
function if there exists a random set ^ (or the corresponding random indi- 
cator function) such that ^ is almost surely measurable and Eg/i(^) = ^{gh) 
for all /i € "^o- By approximating the atomic masses at two points with con- 
tinuous functions, it is easy to see that the weak realisability is equivalent to 
^{Sx,y) = Px,y for 1/ (8) i^-almost all {x, y), in contrast to the strong realisabil- 
ity requiring this equality everywhere. The strong and weak realisability do 
not coincide in general. For instance, a non-positive function which vanishes 
almost everywhere, but takes some negative values is weakly realisable by 
the empty set, but not strongly realisable. Nevertheless, in the case of a 
stationary random regular closed set ^ in M*^ and the Lebesgue measure 
the strong and weak realisability properties coincide, see Theorem 4.7. 

In view of the required continuity property of functions from G it is essen- 
tial to ensure that gh{F), F G T, defined in (4.2) is continuous in the Fell 
topology. Note that it is not the case for most non-trivial measures u, even 
for the Lebesgue measure. The continuity holds only on some subfamilies of 
J-" considered in the following sections. 

4.3. Neighbourhoods of closed sets. For simplicity, in the following con- 
sider random sets in the Euclidean space, i.e. assume that X = R'^ with 
Euclidean metric d. 

Let J-"^ be the family of e-neighbourhoods of closed sets in W^, i.e. J-^ 
consists of = {x : d{x, F) < e} for all F ^ F and also contains the 
empty set. The vector space G is generated by constants and the functions 
gh defined by (4.2) with the Lebesgue measure v. 

Lemma 4.5. The space F^ with the Fell topology is compact and, for 
each h & 'rfo, the function gf^ is continuous on T"^ . 

Proof. Recall that the upper limit of a sequence of sets {-Fn} is the set 
of all limits for sequences such that G Fn^. for all /c, while the 

lower limit is the set of all limits for convergent sequences such that 

Xn € Fn for all n. The sequence of closed sets converges in the Fell topology 
if its upper and lower limits coincide. 

If i^n = Q € converges to F in the Fell topology, then we can assume 
without loss of generality (by passing to subsequences) that F„^o converges 
to Fq, so that F = Fq and F £ F'^. Thus, F'^ is a closed subset of F and so 
is compact, since F is compact itself. 

Consider a non-negative /i € supported by a ball Bji centred at the 
origin with sufficiently large radius R. Fn F in the Fell topology, 
then the upper limit of {F^ n Br) is a subset of {F n Br). Thus, gh{F) = 
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gh{Fn n Bn) < gh{{F n B^Y) for any 5 > and sufficiently large n, so that 
g/j is upper semicontinuous on J-" and so on J-"^. 

In order to prove the lower semicontinuity of gh on assume F^ = 
F^ Q ^ F = F§ and Fnfi Fq. Fix 5 > 0. Then the lower limit of Fnfif^Br+s 
includes Fq n Br. Indeed, if x G Fq n B^, then for x„ € and so 

a^n S for all sufficiently large n. Thus, for sufficiently large n, we have 

(Fo n Br) c (f„,o n BR+if . 

Taking (e — (5)-neighbourhoods of the both sides yields that 

(Fo n BrY"^ c (f„,o n Br+sY c (f„ n . 

If rc G Fq~^ n Br^it+s, then there is a point y G Fq with d(x, y) < e — 5, in 
particular y € Br. Thus, 

Taking r sufficiently large yields that g/i(F„) > gft(FQ~^). Since the interior 
of F equals Us>oFq~^, the Lebesgue theorem yields that gh{FQ~^) — > ghiF) 
as (5 — 7> 0, i.e. gh is lower semicontinuous on J^^ . 

For a non-positive function h with compact support, the result follows by 
applying the above argument to its positive and negative parts. □ 

Theorem 4.6. A function Px,y, x,y G M'^, is weakly realisable by a ran- 
dom closed set ^ with realisations in J^^ for given e > if and only if 

$(gO > inf g;,(F), he%, 
where ^{gh) is given by (4-3). 

Proof. In view of the continuity of g^ established in Lemma 4.5, it 
suffices to refer to Corollary 2.7. □ 

In order to handle random sets with realisations from the space J-^ = 
Ue>o^'^, we need the regularity modulus defined as the infimum of 

e > such that F G F^/" and x(F) = oo if F ^ 7"°. 

Theorem 4.7. For any given r > 0, a function Px,y, a;, y € M.'^, is weakly 
realisable by a random closed set such that Ex(^) < r if and only if 

(4.4) inf MF) - g„(F)] + $(gO < r , he^o, 

where ^{gh) is given by (4-3). If, additionally, px^y is an even continuous 
function of x — y, then px,y is strongly realisable by a stationary random 
closed set ^. 
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Proof. Function x is lower semicontinuous, since {F £ J-" : x(-F) < 
c} = T^/^ is closed for all c > 0. Furthermore, 

(4.5) {FeT: x{F) < g^(F)} C {F e T : x{F) < c} , 

where c = J \h{x, y)\v{dx)v{dy) is a finite upper bound for gh{F). The left- 
hand side of (4.5) is compact, since gh is continuous on T^^^ by Lemma 4.5 
and the right-hand side of (4.5) is compact. Thus, x is a regularity modulus 
and the result follows from Theorem 2.6 and (2.2). 

Note that the regularity modulus x is invariant for the group Q of transla- 
tions of M'^. By Theorem 2.10(ii), ^ can be chosen to be stationary. In order 
to confirm the strong realisability, it remains to show that the covariance 
function of a stationary regular closed random set is continuous. 

Since x(0 is integrable, ^ G J^, so that ^ is almost surely regular closed 
and its boundary 5^ has a.s. vanishing Lebesgue measure. By Fubini's the- 
orem, almost every point x belongs to the boundary of with probability 
zero, and so P{x € d^,} = for all x in view of the stationarity property. 

Let P{x, y £ ^} he the covariance function of ^. Take x,y £ W^, and 
{xn,yn) that converges to (x,y). Since with probability 1, x does not belong 
to 9^, ]Ia;gg is almost surely equal to Ia;gint{5) ^or the interior Int(.^) of ^ and so 
^xnG? almost surely converges to IxGg- Similarly, ly^eg — >■ lygg a.s., whence 
the product converges too Ia;„GC,3/„eC ~^ ^xes.,ye^- The Lebesgue theorem 
yields that P{x„,yn € ^} ^ P{x,y € ^}. Since px^y and P{x,y € ^} are 
both continuous and coincide almost surely, they are equal everywhere. The 
continuity of P{x, y G can be also obtained by referring to a result of [21] 
saying that the capacity functional of each stationary regular closed random 
set is continuous in the Hausdorff metric. □ 

4.4. Convexity restrictions. The family C of convex closed sets in M'^ 
(including the empty set) is closed in the Fell topology and it is easy to see 
that the function given by (4.2) is continuous on C. Corollary 2.7 yields 
that px^y is weakly realisable for a convex random closed set if and only if 

^gh) > mf gh{F) 

for the functional $(gh,) given by (4.3). 

Let V be the convex ring in W^, i.e. the family of finite unions of compact 
convex subsets of W^. For F G P let x{P) be the smallest number k, such 
that F can be represented as the union of k convex compact sets. 

Theorem 4.8. Let $ he linear functional defined by (4-3). Fix any r > 
0. Then there is a random closed set ^ with realisations in V such that 
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Eg/i(0 = ^iSh) for all he% and Ex(0 < r if and only if 
mi [x{F) - gh{F)] + ^gh) < r , he^o- 

Proof. The family Vk of unions of at most k convex compact sets is 
closed in T and so is compact, whence x is lower semicontinuous. It is easily 
seen that gh is continuous on convex compact sets, and so is also continuous 
on Vk- Thus, gh is x-regular and Theorem 2.6 applies. □ 

If X = [0,1], then V is be the family of finite unions of segments in 
[0, 1]. The number of convex components of F C [0, 1] is the variation of its 
indicator function, 

n-l 



X{F) = sup^ I lt,eF - It.+ieF 



i=0 



where the supremum is taken over partitions = to ^ ^ ■ ■ ■ < tn = 1- 
The quantity 

v(F) = sup / ip'{x)dx , 
ifiev^, o<ip<i Jf 

where is the family of differentiable functions on [0,1], captures the 
number of components of F with non-empty interiors, in particular v(F) < 
x{F)- Remark that v is not a regularity modulus, because a set F with small 
y{F) can contain an arbitrarily large number of isolated singletons. 

Theorem 4.9. If Px,y is a function of x,y € [0, 1] such that 



(4.6) sup / px,y if' {x)(p'{y)dxdy = oo , 

then there is no random closed set ^ satisfying Ex(0^ < ^jo having p^^y as 
its two-point covering function. 

Proof. Let H be the family of functions h{x,y) = ip'(x)ip'{y) for ip G 'tf^ 
with < (/3 < 1. Then 

v{Ff = sup / 1x,yeF h{x, y)dxdy = sup gh{F) . 

Theorem 2.6 implies that $ is realisable by a random closed set ^ with 
Ex(0^ < c» if and only if 



sup 



inf[xiFf-gh{F)] + HgH) 



< OO , 
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It implies in particular 



sup 

/igH 



mf[x{Ff-gh{F)] + Hgh] 



< oo . 



Since ^ — &h{F) for h £ H, this condition would imply that 

sup^>(g/j) < oo, 

contradicting (4.6). Thus $ is not realisable. □ 

Further results on realisability of random sets can be found in [5], where 
it is shown that by relaxing the closedness assumption it is possible to split 
the positivity and regularity conditions as it was the case in Section 3.3. 



5. Contact distribution functions for random sets. Results from 
Section 4 concern realisability of the two-point covering probabilities, which 
are closely related to the values of the capacity functional (hitting proba- 
bilities) on two-point sets. Here we consider the realisability problem for a 
capacity functional defined on the family of balls in M*^. If T is the capacity 
functional of a random closed set then 

t{Br{x)) = P{c n Br{x) + 0} 

is closely related to the spherical contact distribution function P{d(x,.^) < 
B\x ^ ^}, -R > 0, which is the cumulative distribution function of the dis- 
tance between ^ and x given that x ^ ^. 

Theorem 5.1. A function Tx{R), R > x £ A C W'-, is realisable as 
T{Bji{x)) for a random closed set ^ if and only if 

m 

(5.1) Hg) = Y.a,T,^{R,)>0 

1=1 

for all m > 1, xi, . . . , Xm € A and . . . , Rm > 0, such that the function 

m 

(5.2) g(F) = ^ a, lB^^(,^)nF^0 > , F G ^ 

1=1 

is non-negative. 
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Proof. The necessity is evident. 

Sufficiency. Let G be the vector space generated by constants and func- 
tions gh^x{F) = h{d{x,F)), F £ T, where d(x,F) is the distance from 
X € to the nearest point of F, and h is a continuous function on R with 
bounded support. The functions gh^x are all continuous in the Fell topol- 
ogy, since the Fell topology in M'^ coincides with the topology of pointwise 
convergence of distance functions d{x,F) for x G M'^, see [20, Th. B.12]. 

It suffices to show that $ is positive on G. Let g{F) = "^^i aihi{d{xi, F)). 
Uniform approximation of hi, . . . ,hm by step functions on their supports 
yields a function g of the form (5.2) so that g{F) > —e for some e > 0. 
Letting e J, and using (5.1) yield that 

m „ 

^g) = Y^ai / hi{t)dTx,{t)>^. 
i=i •' 

□ 

If Tx{R) = t{R) does not depend on x, it may be possible to realise it 
as the contact distribution function of a stationary random closed set. If 
the argument x of Tx{R) takes only a single value, then the necessary and 
sufficient condition on Tx{-) is that it is a non-decreasing right-continuous 
function with values in [0,1], i.e. the cumulative distribution function of a 
sub-probability measure on R+. The following result concerns the case of x 
taking two possible values. 

Theorem 5.2. Let xi,X2 G M'^, with I = \\xi — X2\\, and let Tx^ and Tx2 
be cumulative distribution functions of two sub-probability measures on R_|_. 
Then there exists a random closed set ^ such that Tx^iR) = T{Bji{xi)) for 
r > and i = 1,2 if and only if for all r >0 

(5.3) Tx,imax{R-l,0)) < Tx,{R) < Tx,{R + 1). 

Proof. Necessity: Let ^ be a random closed set with Tx^ {R) = T{Bfi{xi)). 
Let ai and 02 be random points such that oi, 02 € C a-S- and Ri = d{xi, Oj) = 
d(xi,^), i = 1,2, have cumulative distribution functions Tx^ and Tx2 respec- 
tively. Then \Ri — i?2| < Indeed, if, for instance, Ri > R2 + I, then 02 
is nearer to xi than oi contrary to the assumption. Thus Ri < R implies 
R2 < R + I, so that Tx^{R) < Tx^iR + I)- The symmetry argument with xi 
and X2 interchanged yields (5.3). 

SufRciency. Define two random variables Ri and R2 as inverse functions 
to Tx^ and Tx2 applied to a single uniform random variable, so that (5.3) 
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yields that \Ri — -R2I ^ ^ a-S- This means that none of the balls -B/jj(3;i) and 
Bfj^{x2) lies in the interior of the other one. Now construct random closed 
set ^ consisting of two points: ai on the boundary of Brj^{xi) but outside 
of the interior of Br^{x2) and 02 on the boundary of Bii^{x2) but outside 
of the interior of Bfi^{xi). Then ai is nearest to xi and 02 is nearest to X2 
with given distributions of the distance. □ 

APPENDIX: A COMBINATORIAL LEMMA 

Recall that Pt{^) denotes the packing number of X with metric d, i.e. the 
maximum number of points in the space X with pairwise distance exceeding 
t, see [18, p. 78]. 

Lemma A.l. IfY = ^5^^ is a counting measure of total mass n, then 
for all t > 0, 

n 



Proof. Denote 



n{Y,Xi) = Y{Bt{x,))-l, 
where Bt{xi) is the closed ball of radius t centred at Xi. Furthermore, denote 



Then 

g/^.(^-'5xJ=gh,(n-2n(y,x,), 

g;,, (y + = g/,, (y) + 2n(y, X,) + 2 . 

Let Xi and Xj be two distinct points from the support of Y with d[xi,Xi) < t. 
Assume that n{Y, Xi) < n{Y, Xj) or n{Y, Xj) = n{Y, Xj) with i < j and define 

Y' = Y- 6,,^ + 6,^ 

obtained from Y by transferring a mass 1 from Xj to Xj. Call Y" = Y — 5xy 
Remark that n{Y" , Xi) = n(y, Xi) — 1 because d(xj, Xj) < t. Since n(y, Xj) > 
n{Y,Xi), 

= gh,{Y) - 2n{Y,Xj) + 2n(y",x,) + 2 
= gh,{Y) - 2n(y, Xj) + 2n(y, x^) - 2 + 2 
<gh,iY). 
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Furthermore, n{Y', Xj) = n{Y, Xi) because the transferred mass remains in 
the ball with centre Xi and radius t, and n{Y',Xj) = n{Y,Xj) as well. Thus 
n{Y',Xi) < n{Y',Xj). Repeat the mass transfer from xj to Xi until the mass 
at Xj disappears. Call the resulting counting measure Yi. 

Apply the same construction to Yi and repeat it until there are no more 
distinct points at distance at most t. This happens in a finite time because 
the cardinality of the support of Y strictly decreases at each step. 

The obtained counting measure Y is supported by a set of points {yi, . . . , yq} 
with pairwise distances exceeding t. Thus, 

Q 

gh,{Y)>gh,iY) = Y,mi{m,-l), 

i=l 

where = Y{{yi}). Under the restriction rrii = n, the minimal value 
mi{mi — 1) is reached for = n/q, whence 

gh.(n>«(^-l) • 
It remains to note that q < -Pt(X). □ 

It is also possible to define a counting measure by placing masses from 
the interval [n/q, n/q + 1] at the points forming the packing net of X. Thus, 
there exists a counting measure Y such that 
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